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DEPC 
TEC 
TECD 


TFA 
EEA 
SPECO 
AVCO 


Ec 
HRSG 


Nomenclature 

c cost per unit of exergy ($/GJ) 

Ċ cost rate associated with exergy ($/h) 
E exergy rate (kW) 

Edest rate of exergy destruction (kW) 

Eheat rate of exergy transfer by heat (kW) 
f exergoeconomic factor (%) 

h enthalpy (kJ/kg) 

I investment cost ($) 

L Lagrangian 

L rate of the Lagrangian 

IT exergoeconomic cost vector (kW) 

r relative cost difference (%) 

S entropy rate (kW/K) 

Q rate of heat transfer (kW) 

Ww power (kW) 

x vector of decision variables 

y vector of all the inputs and outputs 
Y constraint function 

Z cost rate associated with the sum of capital 


investment and O&M ($/h) 


Abbreviations 
TADEUS thermoeconomic analysis for diagnosis of effi- 


ciency reduction in systems 

diesel engine powered cogeneration 

theory of the exergetic cost 

theory of the exergetic cost—disaggregation meth- 
odology 

thermoeconomical functional analysis 
exergoeconomical analysis 

specific exergetic cost analysis method 

average exergetic cost analysis method 
modified productive structural analysis method 
last in first out principle 

engineering functional analysis 

structural theory of thermoeconomics 
compressor 

air preheater 

turbine 

combustion chamber 

evaporator 

economizer 

heat recovery steam generator 


Greek letters 


n energy efficiency 

E exergy efficiency 

B capital recovery factor (%) 

À vector of all the Lagrange multipliers 
T number of units and junctions 

w number of units, junctions and branching points 
Superscripts 

BQ steam 

CHE chemical 

P mechanical 

T thermal 


W work or electricity 

Subscripts 

av average 

D exergy destruction 

F fuel 

i any stream 

in streams entering a component or system 
k any component 

mar marginal 

out streams leaving a component or system 
P product 

tot total plant 


1. Introduction 


Exergoeconomics (thermoeconomics) combines exergy analy- 
sis with conventional cost analysis in order to assess and improve 
the performance of energy systems [1]. The primary contribution 
of an exergy analysis to the evaluation of an energy system comes 
through a thermoeconomic evaluation which considers not only 
the inefficiencies but also the costs associated with these 
inefficiencies and the investment expenditures required to reduce 
them. The objective of a thermodynamic optimization is to 
minimize the thermodynamic inefficiencies within the system, 
whereas the objective of a thermoeconomic optimization of a 
system is to estimate the cost-optimal structure and the cost- 
optimal values of the thermodynamic inefficiencies in each 
component [2-4]. 

In this paper, starting from earlier pioneer works in the area of 
thermoeconomics, major contributions are summarized. By 
detailed search with appropriate keywords, a researcher can 
easily reach to the related papers in the open literature. First, a 
short historical overview on the exergoeconomic analysis and 
optimization of thermal systems is given. In the following sections, 
exergy costing and cost accounting concepts are described. Various 
thermoeconomic methodologies are discussed and exergoeco- 
nomic analysis and optimization of combined heat and power 
production are presented. 


2. Historical overview 


The idea of coupling exergy and cost streams was first 
discussed by Keenan in 1932 [5]. He pointed out that the value 
of the steam and the electricity rests in the “availability” not in 
their energy [6-8]. In the late 1950s, the studies of second law 
costing started in two different places independently. Tribus and 
Evans [6] studied desalination processes by exergy analysis, 
which led them to the idea of exergy costing and its applications 
to engineering economics, for which they coined the word 
“thermoeconomics” [9]. The concept of their procedure was to 
trace the flow of money, fuel cost and operation and amortized 
capital cost through a plant, associating the utility of each 
stream with its exergy. In the late 1960s, El-Sayed connected 
with Evans and Tribus in their research on desalination and they 
published one of the important papers in the subject in 1970, in 
which the mathematical foundation for thermal system opti- 
mization was given [10]. Another study on their idea was 
performed by Reistad who applied the method of El-Sayed and 
Evans to a simple power plant, comparing that approach with 
conventional optimization procedures [11]. Also in 1960s, 
Gaggioli studied the optimal design of power plant steam 
piping and its insulation in his Ph.D. thesis. He proposed costing 
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steam exergy at a value to that of power produced, penalizing 
exergy destructions and losses for the electricity which there- 
fore would not be produced [12]. 

In Europe, many important works on the second law analysis 
methodologies and on “exergy” itself were performed in late 1950s 
and during 1960s and 1970s [13-20]. Bergmann and Schmidt 
assigned costs to the exergy destruction in each component of a 
steam power plant and optimized feed water heaters [17]. Szargut 
used exergy costing procedure in the analysis of a simple 
cogeneration plant by derivation of an approximate generalized 
formula and introduced ecological cost coefficient into the 
literature. The use of this coefficient made the determination of 
cumulative consumption of non-renewable natural resources in a 
production process possible [18-20]. Some of these works were on 
thermoeconomic analysis and collected in a book in 1985 by Kotas 
[21], which is still accepted as one of the basic references in exergy 
analysis and thermoeconomics of thermal systems. 

Although there had been numerous studies and theoretical 
approaches in the field of thermoeconomics, the methodological 
and functional applications of it to the analysis, design and 
optimization of thermal systems comprehensively did not start 
until 1980s. During this decade, the interest and studies related to 
exergoeconomic analysis techniques and applications have 
increased considerably. Frangopoulos [22-27] and von Spakovsky 
[26-30] applied and formalized the first autonomous method 
developed by Evans and El-Sayed [10]. Tsatsaronis introduced 
some key concepts of thermoeconomics such as Fuel and Product 
[31]. Valero et al. published another key paper on thermoeco- 
nomics and presented the basic methodology related to exergy 
based cost analysis and applications [32] and with Lozano, he 
presented basics and several applications of the theory of exergetic 
cost, a major cornerstone approach to the field of thermoeco- 
nomics [33]. 

The major contributions were done in 1990s, to achieve a 
greater standardization and formalism in the area of thermo- 
economic studies. The common idea was to propose a standard and 
common mathematical formulation for all thermoeconomic 
methodologies employing thermoeconomic models that can be 
expressed by linear equations [34]. At that time, Tsatsaronis 
proposed the term “exergoeconomics” which was defined as a part 
of thermoeconomics. Since the latter have been used in a general 
sense expressing the interaction between any thermodynamic 
variables and economics, he suggested that exergy based cost 
accounting methodologies should be indicated by exergoeco- 
nomics [35]. One of the most interesting works was the CGAM 
problem [36-40] developed originally by Frangopoulos, Tsatsar- 
onis, Valero and von Spakovsky. The CGAM problem was named 
after the first initials of the participating investigators. The 
objective of the problem was to show how the methodologies 
were applied, what concepts were used and which results were 
obtained in a simple and specific problem [36]. The sequential 
papers published in the same issue had the common definitions of 
the physical, thermodynamic and cost models plus the objective 
function but different procedures of thermoeconomic optimiza- 
tion. The authors emphasized that the aim of the presented work 
was the unification of thermoeconomic methodologies, not a 
competition among them. Another important project called 
TADEUS (in honor of Tadeus Kotas [21]) was initiated in 2001 
[41-43]. The major goals of that project were to apply procedures 
from different research groups in thermoeconomic analysis to the 
diagnosis of the energy system malfunction and inefficiencies and 
to establish the common concepts and nomenclature and compare 
the results and highlight the main characteristics of each approach. 

Since 1980s, there have been numerous published papers all 
around the world on exergoeconomic cost analysis, application 
and optimization of thermal systems. Good fractions of them have 


been published since mid-90s due to the improved structural 
formalism of the exergoeconomic methodologies [44-81]. Exer- 
goeconomic studies have been successfully applied to power 
plants, combined heat and power production or cogeneration 
facilities [45-47,50,51,54,57,59-81]. Progressed innovative meth- 
odologies of fuzzy logic and genetic algorithm methods have also 
been applied to existing power plants and cogeneration facilities 
since late 1990s [47,59,69]. 


3. The cost of exergy and exergoeconomics 


The exergy cost of a mass or energy stream is the amount of 
exergy required to produce it. In the case of a cogeneration plant, 
the exergy cost of the net power is the exergy provided by fuel to 
generate the electrical power delivered to the network by the 
cogeneration plant [33-39]. When analyzing the cost structure of a 
system, one can distinguish between average costs and marginal 
costs [1-3,33]. Average costs are ratios and express the average 
amount of resources per unit of product, and can be defined 
mathematically as [2,8,21,33] 


Cavi 


É, (1) 


Cavi = 
where Cav,i is the unit specific average cost of the ith stream, Gus 
and E; represent the exergetic cost and exergy of the ith stream, 
respectively. Marginal costs are derivation and indicate the 
additional resources required to generate one or more units of 
product under specified conditions. Mathematically, it is defined as 
[33] 


ac; 


aE, (2) 


Cmari = 

Since the average cost is not predictive, it can only be known 
after production process, when we know how much resources 
were used. On the other hand, marginal costs are predictive in 
nature; they can be used to calculate additional fuel consumption 
when the production is modified. Thermoeconomic optimization 
methods are generally based on the marginal costs while average 
costs are used in the exergoeconomic analysis of the systems [32- 
34]. In exergy costing, a cost is associated with each exergy stream. 
Exergy transfers by the entering and exiting streams of matter and 
by power and heat transfer rates may be written, respectively as 
[2] 


Ġ = cÈ; (3) 
Co = CeEe (4) 
Cw = CwW (5) 
Cy = CqEq (6) 


Accordingly, for a kth component receiving heat transfer and 
generating power, we may write [2] 


XO (CeEe)g + Cw Wk = Cq pEge + >> (CE), + (7) 


e i 


In the case of a single product of the plant, the unit cost of the 
product can be easily determined using Eq. (7). In the case of a 
multi-product plant such as cogeneration plants, the cost-balance 
equation is not sufficient. Additional criteria are required to 
determine the relationship between the unit costs of the different 
products. This is where exergy can be used as a basis for cost 
allocation of the products. 

A diesel engine powered cogeneration (DEPC) plant (see 
Fig. 1) is considered as an illustrative example of the exergy 
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— Fuec] oil 


Water 


Fig. 1. The Schematic of an actual diesel engine powered cogeneration plant. C: compressor, T: turbine, WHB: waste heat boiler, DeSO,: desulphurization, AWR: air-water 
radiator, IC: intercooler, LOC: lubrication oil cooler, LOT: lubrication oil tank, SD: steam drum, FWT: feed water tank, FDT: fuel oil day tank, CON: condenser, FFM: fuel 


forwarding module, and P: pump. 


costing of co-generated products. When the engine starts, air is 
charged to the compressor of the turbocharger unit. The 
turbocharger consists of a turbine and a compressor, which 
are mounted on a common shaft. The shaft power obtained by 
the operation of diesel engine is transferred to the generator for 
electricity production. The exhaust gases leaving the engine flow 
through the turbine of the turbocharger unit to produce the 
necessary shaft work for the compressor. The air leaving the 
compressor is cooled by water in an intercooler before air enters 
the engine cylinders. The exhaust gases leaving the turbine enter 
the exhaust gas boiler unit to transfer heat to the feed water to 
produce steam for manufacturing facilities in the factory and for 
preheating of streams in the auxiliary equipments such as fuel 
forwarding module (FFM) and fuel oil in daily usage tank (FDT). 
The exhaust gas leaving the boiler is sent to the DeSO, 
(desulphurization) unit in which the SO, emission is lowered 
to the acceptable legal values. After the DeSO, unit, the exhaust 
gases are released to the atmosphere. In Table 1, exergetic cost 
rate balances and related auxiliary equations for each subsystem 
of the DEPC plant are presented [80,81]. 

The explicit determinations of auxiliary equations based on 
exergy criteria were presented in previous works [7,8,12]. 
These methods are known as the equality, the extraction and 
the by-product methods. All these methods require a 
judgment regarding the purpose of each unit of the plant and 
distribute the costs proportionally to exergy contents of the 
flow. These three methods are applied to the turbine of the 
turbocharger unit of the DEPC plant (see Fig. 1) as an illustrative 
example. The exergetic cost-balance equation of the turbine 
from Table 1 is 


Ce — Cy + Zrursine = C Weine ee 
or 
cgE¢ + Zrursine = C7E7 + CwEw (8b) 


Assuming the exergetic cost of the inlet stream ce is known, the 
cost-balance equation has two unknowns, namely, c7 and cw. Thus, 


an assumption must be made regarding the allocation of the cost to 
the two outputs: shaft power Ew and the exit stream E7. 


(i) In the equality method, the generation of the two products is 
considered to have the same priority, so the cost of high 
temperature exhaust gas of the turbine is charged to the 
production of both outputs (products) proportionally to their 
exergies as 


c =w= C6E6 + Żrursme (9) 
E7 + Ew 
(ii) In the extraction method, it is considered that the purpose of the 
turbine is to generate shaft power and thus the entire cost of the 
turbine must be charged against it. This results in assuming the 
unit exergetic cost of the exhaust gas entering the turbine to be 
the same as that of the exhaust gas leaving it. Thus, 


C7 = C6 (10) 


_ Ce(És — £7) + Zrursine (11) 


Cw 
Ew 


(iii) In the by-product method, the cost of one of the outputs is 
assumed to be known. It is assumed that the generation of 
process stream is essential even if no shaft power would be 
generated. Therefore, the exergy of the exhaust gas is costed as 
if it were produced alone in a diesel engine at the required 
pressure and temperature. Therefore, the cost of shaft power is 
determined by the cost balance. 


Table 1 shows the exergetic cost rate balances and correspond- 
ing auxiliary equations for each subsystem of the DEPC plant given 
in Fig. 1 [80,81] using the extraction method. 


3.1. Exergetic cost formation process and cost accounting 
methodologies 


3.1.1, Fuel-product model 
Exergetic cost accounting provides a wide and clear vision of the 
usage and degradation of energy and in consequence of natural 
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Exergetic cost rate balances and corresponding auxiliary equations for each subsystem of the DEPC plant. State numbers refer to Fig. 1. 


Component 


Exergetic cost rate balance equations for the components of DEPC 


Auxiliary equations 


Compressor Cites + Zcomp = Ġ — G 


Intercooler (Coo — C1) + (G3 — Goa) + Zc 


Ġ =0 (É =0) 
Cy = Go aes Ca = G3 


G- : a ae : 
E21 — E29 E24 — E23 


Lubrication oil cooler (Coa — Cos) + Zroc = Cog — C27 Eng pz E 
C24 = C25 
Go _ Cr 
Diesel engine C3 + Cs + (C21 — G2) + (C26 — G7) + Zoe = Co + Cw-eectric eee att 
G5 _ Gy 
E26 E27 
Turbine Co — Ċ7 + Zrursine = CWrypeine Ce =C7 
eee ae Ga _ G3 
Waste heat boiler (G — Cs) + Zwe = G3 — Ci2 Ey Ég 
C7 =Ce 
DeSO,, unit (Cy — Cg) + Zpeso, = 0 Co =0 
F : z E A Gs _ Ge 
Fuel oil day tank (Gis — C16) + Zepr = C4 — Coro === 
Eis E16 
Fuel forwarding module (G7 — Gs) + Êm = Cs — Cy Cu _ rp _ Gis 
Ei7 Fig Eie 
Nee Gs _ Co 
Condenser (C9 — Go) + Zcon = Cig — Cis lim, Be 
C299 = 0 
5 bes 7 Go € 
Pump Wome + Ze = Cu — Go zal) = Sal 
Eio En 
ONR : A f ee: Cu Ce © 
DEPC Ci + C5 + Go + Zperc = C14 + Cis + Co + Cw-ELECTRIC a4 ob 2 
Eig Eis Er 


resources [34]. Cost accounting methodologies are based on cost 
allocation rules. In order to identify the process of cost formation, 
general relationships are required that relate the overall efficiency of 
the plant and the cost of the products with the efficiency and the 
irreversibilities of each individual component that form the system 
[33,34,48,61 J. 

In order to build up the cost formation in a system, productive 
structure information of the system must be obtained. This 
structure consists of a qualitative and quantitative analysis of the 
relations that link the flows of exergy inputs and outputs between 
components of an economic unit. Since all thermal plants have a 
defined aim to produce one or more products, the quantity of 
resources must also be identified through mass or energy flows 
which are known as fuel. Thus, each of the components of the plant 
can be characterized by its fuel and its product(s). In Table 2, the 


Table 2 
Definitions of the exergies of the fuels Ép and the exergies of products Ep for the 
components of the DEPC plant. State numbers refer to Fig. 1 [82]. 


Component Ér Ep 

Compressor Ec E Èi 

Intercooler EE (Eo1 — Exo) + (Eoa — Eo3) 
Lubrication oil cooler jog ER Boys = Dyn 

Diesel engine ES Aalia delish Bb os Bop oh Bg 
Turbine Eg — É Er 

Waste heat boiler É —Eg EEIE 

Fuel oil day tank EREE Bin = Bare 

Fuel forwarding module EPE [Bis ESE 

Condenser ERE EIS EEE 

Pump Ep E11 — E10 


fuel and product exergy definitions of each component in the DEPC 
plant is given [82]. 

In accordance with the fuel-product model, the product of one 
component is used as fuel for another component or as part of the 
total production of the plant. Thus, 


n 
Epi = Erio + > Érij i=0,1,...,n (12) 


jal 


where Erij is the production portion of the ith component and is 
the fuel of the jth component. In Eq. (12), the component “0” is 
considered as the environment, and then E;;9 represents the 
production portion of the ith component which leads to the final 
product, coming from the environment to the ith component. The 
resources entering each component can also be expressed as 


n 
Eri =Eroit+ > Érji i=0,1,...,n 
= 


(13) 


where Ér o.i represents the external resources entering to the plant, 
which enters through ith component. The total fuel and product of 
the system can be obtained as [34-36,39] 


n 
Er total = > ÉF oj (14) 
jal 
. n . 
Ep total a X Ér jo (15) 
j=l 


3.1.2. The cost equations of fuel-product model 
The proposed fuel-product model can be applied to exergetic 
cost structure formation of a system and exergetic costs of product 
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and fuel can be expressed, respectively as 


Ĉi = Cio + > Gray i=0,1,...,n (16) 
jal 

. . n . 

ri = Gort $ Čj i=0,1,...,n (17) 
J 


The proposed cost model of fuel-product approximation is the 
base of the “Exergetic Cost Theory” [33]. It is explained in the next 
section with the cost assessment rules of the theory. The exergetic 
cost-balance equations and auxiliary relations of the DEPC plant 
given in Table 1 are based on the fuel-product model and obtained 
by using SPECO approach that is one of the important thermo- 
economic methodologies described in the next section. 


4. Thermoeconomic methodologies 


Following El-Sayed and Gaggioli [7,8], thermoeconomic meth- 
ods can be classified in two groups: algebraic methods and calculus 
methods. These two basic categories are examined with their 
subcategories and related works. 


4.1. Algebraic methods 


The algebraic methods use algebraic cost-balance equations 
derived from conventional economic analysis and auxiliary cost 
equations for each subcomponent of any system presented [2]. 
They are related with the cost formation process of the system in 
order to investigate the average costs. 


4.1.1. The theory of the exergetic cost (TEC) 

Theory was developed by Lozano and Valero [33] and the 
methodology presented in theory is based on a set of propositions 
[2]. This theory permits the introduction of a new thermodynamic 
concept called exergy cost. For given a system whose limits, 
disaggregation level, and production aim of the subsystems have 
been defined, exergy cost E* of a physical flow is defined as the 
amount of exergy needed to produce this flow [33]. The first step in 
the application of this methodology is the division of the system into 
units, which may adapt to a component or a set of components. A 
single product and fuel for each component in the system must be 
defined. Thus, a system of equations can be built with a cost-balance 
equation for each unit (proposition 1), cost allocation equations for 
external flows into the system for which costs are externally defined 
(proposition 2), and losses for which the cost is set equal to zero 
(proposition 3). However, this procedure is rarely enough. The 
solution may be possible when two following propositions are 
considered: (1) if definition of fuel of acomponent includes a stream 
that goes through another component and is used in it, then the unit 
cost of the stream flowing into and out of the component is the same 
(proposition 4a). (2) If the product of a component is composed of 
two or more streams then the unit cost of those streams are equal 
(proposition 4b). Using the F (fuel)-P (product)-L (loss) definitions, 
corresponding matrices can be developed. Starting from these 
matrices and using the data from design and operation of the any 
thermal power plant, it is possible to carry out the exergetic and 
energetic analysis of the plant. The unit exergetic costs obtained 
from the procedure can be used for the optimization of any thermal 
plant. This methodology is straightforward in its application, 
allowing for improvements and modifications. 


4.1.2. The theory of exergetic cost—disaggregating methodology 
(TECD) 

This theory was proposed by Valero and Lozano [33,34] 
searching the thermoeconomic unification methodology which 


is similar to Frangopoulos’ Thermoeconomical Functional Analysis 
(TFA) [22-25]. The exergetic cost determined is distributed among 
all other components according to the variation of entropy within 
each one. 


4.1.3. Exergoeconomic analysis (EEA) methods 

The exergoeconomics approach was proposed by Tsatsaronis 
and co-workers [1-4,35-37,52]. There are two possible variants in 
these methodologies: specific cost and average cost. The average 
cost concept is very similar to the application of Theory of the 
Exergetic Cost. In the specific-cost method, the cost of the addition of 
exergy to a material/stream current is determined and charged to 
the unit that makes use of that exergy. This means that a 
component will obtain the exergy from a stream at different costs, 
depending on the components that supplied the exergy to that 
stream. Due to the purpose of the energy conversion plant (i.e. 
power plant or combined heat and production plant) the division 
between subcomponents of the main plant may be made. If power 
produced is the primary product, cost of all external irreversi- 
bilities is loaded on the electricity. If heat is produced as primary 
one, the similar procedure is followed by the steam produced in the 
plant. The optimization procedure of the exergoeconomics is based 
on an iterative design improvement procedure that does not aim at 
calculating the global optimum of a predefined objective function 
[35], as the conventional optimization methods perform, but tries 
to find a good solution for the overall system design. The methods 
involve some characteristic parameters of exergoeconomic ana- 
lysis such as relative cost difference, exergoeconomic factor, and 
exergetic efficiency. These methods can be grouped in three 
classes: (i) Last-in-First-out (LIFO) principle, (ii) Specific Exergy 
Costing/Average Cost (SPECO/AVCO) approach, and (iii) Modified 
Productive Structure Analysis (MOPSA) approach. 

Last-in-First-out (LIFO) principle. This approach was developed 
by Tsatsaronis et al. [1]. LIFO principle of cost accounting is used for 
the spent exergy units to calculate the cost of exergy supply to the 
carrier. The method eliminates the need for auxiliary assumptions 
in the exergoeconomic analysis of energy systems and improves 
the fairness of the costing process. This approach is based on a 
simple idea, which is directly linked between asset management 
and valuation method: The exergy units that are supplied last to 
the material stream are used first. Thus, it is straightforward to 
calculate the cost associated with the exergy units removed from a 
material stream in a process step provided that this cost is 
calculated in the previous steps during which the exergy units 
currently removed from the stream are then supplied to it. 

Specific Exergy Costing/Average Cost (SPECO/AVCO) Approach. 
This method was first presented by Lazzaretto and Tsatsaronis 
[48,52]. Following improvements of the methodology by Tsatsar- 
onis and co-workers it includes combination with fuzzy inference 
systems for a more exact exergoeconomic evaluation of plant 
components [54] and systematic improvement of SPECO metho- 
dology [71]. According to the approach, fuels and products are 
defined systematically by registering exergy additions to and 
removals from each material and energy stream. The records of 
cost additions to and removals from the same stream in 
conjunction with the application of the LIFO principle are then 
used to calculate local average costs. This method consists of three 
main steps: (1) identification of exergy streams, (2) definition of 
fuel and product for each system component, and (3) cost-balance 
equations. This method has been largely and successfully used and 
applied to thermal systems by the researchers in the area of 
thermoeconomics [4,67,71-81]. 

Modified Productive Structure Analysis (MOPSA) approach. This 
approach was first presented by Kim et al. [45]. The method is 
based on the exergy costing approach without flow-stream cost 
calculations. Lozano and Valero [33] presented the cost evaluation 
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technique of the thermal system, which provided a theoretical 
basis for the exergy costing method. In the MOPSA method, a cost- 
balance equation is obtained by assigning a unit exergy cost to each 
disaggregated exergy in the stream at any state. The monetary 
evaluations of various exergy costs, as well as the production cost 
of electricity for the thermal system, are obtained by solving the set 
of equations for the unit of exergoeconomic costs. This method has 
been applied and compared with SPECO approach for the 
exergoeconomic analysis and optimization of combined heat 
and power plants [45-57,60]. 


4.2. Calculus methods 


Calculus methods are built on differential equations. Cost flows 
in a system are developed in a link between optimization 
procedures that are based on the Lagrange multipliers and they 
are used to determine marginal costs. The characteristic difficulty 
in the application of calculus methods to complex systems such as 
combined heat and power plants is the fact that the Lagrange 
multipliers vary from iteration to iteration when component 
thermoeconomic isolation is not achieved [6]. Following 
approaches are used in calculus methods. 


4.2.1. Thermoeconomic functional approach (TFA) 

This approach was first studied as a Ph.D. thesis by Frangopoulos 
[22-27] while the first remarkable application of it was the CGAM 
problem [38]. This methodology is based on the Lagrangian method 
of mathematical optimization. Its complete implementation 
requires the existence of a sufficiently accurate simulation of the 
system in order to determine the first order derivatives of the 
objective function. The method is based on the decomposition of the 
system into subcomponents, which may or may not correspond toa 
physical component of the system. Each component has a single 
product, and costs can be determined by the solution of the cost 
balances of the system components. The TFA optimization is based 
on the direct use of an (usually nonlinear) optimization algorithm 
since it requires the least effort in the analysis of complex systems. 
The method introduces the functional diagram concept and assigns 
only one function and one product to each component so that 
auxiliary equations are not needed in the solution procedure. This 
method does not give any information about the physical and 
economic interrelationships among the system components. On the 
other hand, scaling of the variables and of the objective function is 
required to achieve convergence to the optimum point. Discrepan- 
cies that exist in searching convergence to the global optimum can 
be solved by searching different initial points in the problem. In 
order to find the optimum solution of more complex systems, 
Frangopoulos formulated Intelligent Functional Approach (IFA) which 
was a further development of the TFA [24,25]. 


4.2.2. Engineering functional analysis (EFA) 

The basics of the EFA theory were developed by von Spakovsky 
and Evans [30]. The first laborious work of the method was the 
CGAM problem [39]. The basic decomposition algorithm of the EFA 
method is the variation on the algorithm proposed by Frangopoulos 
[22,23]. It is based on a multi-dimensional version of Modified 
Regula-Falsi Method [83]. In the EFA approach, any thermoeconomic 
model exists on two levels: as a system basis model and as a set of 
detailed subgroup models [39]. These two models include informa- 
tion on the internal geometry and material composition of each 
subgroup. The extent to which these individual subgroup optimums 
are consistent with the global system optimum depends on how well 
the isolation of the subgroups is established by the optimization of 
the system basis model. The optimization procedures of two 
interrelated models mentioned is possible with a certain number of 
iterations which occur in the link between the models and internal 


economy of the system presented. The decomposition strategy used 
in calculus methods in thermoeconomics is based on the Principle of 
Thermoeconomic Isolation by Evans [6]. The idea is stated as follows: 
A component of a thermal system is thermoeconomically isolated 
from the rest of the system if its production P; and the unit cost of the 
resources À; are known quantities and independent from the rest of 
the component variables. Decomposition may only approach the 
global optimum because of the degree of thermoeconomic isolation 
of the independent variables, the choice of subgroups and their 
functions, and the nature of dependent variables [83]. 


4.2.3. The structural theory of thermoeconomics (STT) 

This was proposed as a standard and common mathematical 
formulation for all thermoeconomic methodologies employing 
thermoeconomic models that can be expressed by linear equations 
[34]. The Theory of Exergy Cost (TEC), the SPECO/AVCO approach 
and the Thermoeconomic Functional Analysis (TFA) can be dealt 
with the structural theory. The Last-in-First-Out (LIFO) approach 
can also be reproduced with the structural theory. 


5. Comparison of the thermoeconomic methodologies on the 
CGAM problem 


The CGAM problem was constructed to examine the optimum 
design parameters of a simple combined heat and power plant 
which produces 30 MW electricity and 14 kg/s of saturated steam 
at 20 bars [36]. The two optimization criteria were the maximiza- 
tion of exergetic efficiency and the minimization of the total cost 
rate of the plant. The installation of the plant consists of a gas 
turbine followed by an air preheater which uses part of thermal 
energy of the gases leaving the turbine, and a heat-recovery steam 
generator in which the required steam is produced (see Fig. 2) 
[2,36,53]. The base case design data of CGAM problem is given in 
Table 3 [2,36-39,53]. 

In the definition of CGAM problem, three models were 
considered: the physical model which describes the behavior of 
the system; the thermodynamic model which describes the 
equations of state used to calculate the thermodynamic properties; 
and the economic model which describes the equations for 
calculating the capital costs of the components. In order to simplify 
these models, the following assumptions were made: (i) the air and 
the combustion gases behave as ideal gases with constant specific 
heats, (ii) for combustion calculations, the fuel is taken to be 
methane CHy, and (iii) all components except the combustion 
chamber are adiabatic. 


5.1. Conventional approach 


In the first study of the CGAM problem [36], physical, 
thermodynamic and economic models were described in detail 
for each component of the plant (see Fig. 2). In this primary study, 
the optimization problem was expressed as the minimization of 
the total cost flow rate of the plant Cr, subject to the constraints 
imposed by the physical, thermodynamic, and cost models of the 
installation. That is, 


Cr = cpt FLHV + Zc + Lap + Zec + Zr + Zursc (18) 


where cris the cost per unit energy of the fuel for the total plant, m p 
is the mass flow rate of the fuel, LHV is the lower heating value of 
the fuel, and Z, is the total capital cost associated with the 
subcomponents of the plant, i.e., compressor, air preheater, 
combustion chamber, turbine and heat recovery steam generator. 
This first solution of the CGAM problem was obtained through 
conventional optimization techniques. The decision variables of 
the CGAM problem were the compressor pressure ratio P>/P,, 
isentropic compressor efficiency nsc, isentropic turbine efficiency 
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Process steam 


Feed water 


Exhaust gas 


Net Power 


Fig. 2. The flow schematic of the CGAM problem. C: compressor, AP: air preheater, CC: combustion chamber, T: turbine, Ev: evaporator, Ec: economizer, ST: steam tank, HRSG: 


heat recovery steam generator. 


Ns» temperature of the air entering the combustion chamber T3, 
and temperature of the combustion products entering the turbine 
T4. All other thermodynamic variables were calculated as functions 
of the decision variables. The mathematically optimal solution 
results of this introductory study [36] are given in Table 4. 


5.2. Exergoeconomic methodology approach 


In the second study of the sequential CGAM problem papers 
[37], exergoeconomic techniques were applied and discussed for 
design optimization. The exergoeconomic evaluation of the system 
considered (see Fig. 2) consists of a detailed exergy analysis, an 
economic analysis conducted at the component level, calculation 
of the cost of each stream using an appropriate exergy costing 
method [31], and a final evaluation at the component level 


Table 3 
Base case design data of the CGAM problem [2,36-39,53]. 


State Substance Mass flow rate (kg/s) Temperature (K) Pressure (bar) 


1 Air 91.2757 298.150 1.013 
2 Air 91.2757 603.738 10.130 
3 Air 912757 850.000 9.623 
4 Exhaust gas 92.9176 1520.00 9.142 
5 Exhaust gas 92.9176 1006.162 1.099 
6 Exhaust gas 92.9176 779.784 1.066 
7 Exhaust gas 92.9176 426.897 1.013 
8 Water 14.0000 298.150 20.000 
9 Water 14.0000 485.570 20.000 
10 Fuel (methane) 1.6419 298.150 12.000 


Table 4 
Optimum values of the decision variables and costs in the CGAM 
problem through mathematical optimization method [36]. 


Decision variables Optimum values 


P2/P, 8.5234 

Nec 0.8468 

Nst 0.8786 

Tə (K) 914.28 

Ta (K) 1492.63 
Total cost rate ($/s) 0.362009 
Fuel cost rate ($/s) 0.325489 
Investment cost rate ($/s) 0.036520 
Cost of combustion chamber ($) 0.1469 x 10° 
Cost of air compressor ($) 0.1348 x 10’ 
Cost of gas turbine ($) 0.1927 x 107 
Cost of air preheater ($) 0.8277 x 10° 
Cost of HRSG ($) 0.1202 x 107 


conducted with the aid of some relevant exergoeconomic variables 
[2]. The exergy analysis was performed to identify the sources of 
the thermodynamic inefficiencies and to facilitate development of 
design changes to improve the overall system efficiency. Thus, the 
exergy destruction in a component can be calculated from the 
exergy balance as 


: <TOT <TOT 
Ep = X Éim 3 X Éon (19) 
where E is the flow rate of total exergy associated with the ith 
material or energy stream at the inlet (for subscript in) and outlet 
(for subscript out) of the component. The exergy values of fuel and 


product for the kth plant component, exergy destruction may be 
written as 


Épp = Ere — Épk — Exp (20) 


where E,, represents the exergy loss in the kth component. In 
these papers, the performance of an energy system from the 
thermodynamic point of view was evaluated in addition to the 
exergy destruction, the exergetic efficiency, and two exergy 
destruction ratios [2,31,37]. The exergetic efficiency is the ratio 
of product exergy to the fuel exergy: 


Ek = Erk (21) 
Ere 
In addition to the exergy destruction Ep, and the exergetic 
efficiencies €k, the exergy destruction ratio ypx is used in the 
thermodynamic evaluation of a component. This ratio compares 
the exergy destruction in the kth component with the total fuel 
exergy supplied Ér totai to the overall system: 


È 
Ws (22) 
F total 


Alternatively, the exergy destruction rate of the kth component 
can be compared to the total exergy destruction rate Ep total: 


E 
You = ep (23) 
Ep total 
The exergy loss ratio is defined similarly to Eq. (22), by 
comparing the exergy loss to the total fuel exergy supplied to the 
overall system 


E 
YL.total = —- (24) 
EF total 
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Table 5 


Comparison of the global optimum values of decision variables and costs with the two workable design condition optimal values in the CGAM problem through 


exergoeconomic analysis and optimization methodology [37]. 


Decision variables Optimum case Case 14 [37] Case 22 [37] 

P>/P, 8.52 10.0 8.0 

Nse 0.847 0.850 0.850 

Nst 0.879 0.875 0.870 

T3 (K) 914.3 900.0 933.0 

Ta (K) 1492.6 1485.0 1500.0 
Component Ar (%) Ae (%) Ar (%) Ae (%) Ar (%) Ae (%) 
Air compressor 4.08 1.30 9.04 1.82 4.18 1.31 
Air preheater 22.03 —5.64 12.73 —4.27 25.28 —5.93 
Gas turbine 4.56 0.87 1.35 0.50 2.88 0.73 
HRSG 36.71 —15.26 38.07 —15.58 37.02 —15.30 
Dy pat ($/h) 46.70 47.4 49.4 

Di tot ($/h) 26.70 26.8 27.2 
Estimated value of the objective function ($/h) 1302.8 1307.9 1307.6 
Calculated value of the objective function ($/h) 1303.4 1308.8 1305.3 


The detailed economic analysis of the CGAM problem presented 
in the second study [37] was conducted at the component level as 
part of the exergoeconomic evaluation and simplified by assuming 
some economic factors such as capital recovery factor (6), annual 
system operation time at the nominal capacity (t) to specific 
numerical values. Since the exact formulations of the cost balances 
and additional auxiliary equations depend on the exergy costing 
method applied (see Eqs. (1) through (7)) and on the types of 
exergy values (i.e., total, chemical, physical, thermal, and 
mechanical exergy), two separate exergy costing methods were 
applied [37]. These methods were the traditional average-cost 
exergy costing method (AVCO) [7,8] and the specific-cost exergy 
costing method (SPECO) [48,52]. AVCO method was the simplest 
method applied to CGAM problem among the others. The SPECO 
method was based on a detailed accounting of the cost at which 
every exergy addition to a material stream occurs. Thus in a 
process during which exergy was removed from a stream, the cost 
at which the removed exergy had been previously supplied to the 
stream was calculated and was charged to the product of the 
process [48,52,54,71]. These two methods were applied to the 
problem combined with (i) physical and chemical exergy values 
and (ii) thermal, mechanical and chemical exergy values. In this 
paper, two major exergoeconomic variables were defined for the 
first time: the relative cost difference factor (Ar) and the 
exergoeconomic factor (f) which would be very useful parameters 
in the evaluation of exergoeconomic analysis and optimization of 
thermal systems later [37,45,46,51,52,54,60,73,77,80,81]. The 
relative cost difference rx for the kth component is defined as 


Cpr —C 
Tk = P,k F,k (25) 
Crk 


In an iterative cost optimization of a system, if the cost of fuel of 
a major component changes from one iteration to the next, the 
objective of the cost optimization of the component should be to 
minimize the relative cost difference instead of minimizing the 
cost per exergy unit of the product with this component [2]. In 
evaluating the performance of a component, we want to know the 
relative significance of each category. This is provided by the 
exergoeconomic factor fg defined for the kth component as 


Êy 


fk = = = = 
Zk + Cre(Epk + Erk) 


(26) 


The total cost rate causing the increase in the unit cost from fuel 
to product is given by the denominator in Eq. (26). Accordingly, the 
exergoeconomic factor expresses as a ratio the contribution of the 
non-exergy related cost to total cost increase. A low value of the 
exergoeconomic factor calculated for a major component suggests 


that cost savings in the entire system might be achieved by 
improving the component efficiency (reducing the exergy destruc- 
tion) even if the capital investment for this component will 
increase. On the other hand, a high value of this factor suggests a 
decrease in the investment costs of this component at the expense 
of its exergetic efficiency [2]. Three optimal cases were compared: 
global optimum case of the CGAM problem and two different initial 
workable design case optimal solutions, namely cases 14 and 22 
[37]. Comparison of the results for these three cases is given in 
Table 5. The exergoeconomic analysis and optimization technique 
developed and applied to CGAM problem in the second study of the 
sequential papers show that the method is a valuable tool in the 
exergoeconomic analysis of the complex existing thermal systems 
and it has been applied to different thermal systems in the 
numerous papers [1,4,31,33,37,45,46,51,52,54,60,73,77,80,81] 
since it was published. 


5.3. Thermoeconomical functional analysis approach (TFA) 


The application of the thermoeconomical functional approach 
(TFA) to the CGAM problem was first studied by Frangopoulos [38]. 
He applied three different analysis and optimization procedures in 
his study: (i) direct use of a nonlinear programming algorithm, (ii) 
thermoeconomic functional approach, and (iii) modular simula- 
tion and optimization of the system. In the direct application of the 
nonlinear programming method, the objective of the optimization 
was defined as the minimization of the total cost rate (i.e. capital 
and operation expenses) of the system 


n 
min F = S°Z, + crit LHV, (27) 
z 


r= 


where Z, represents amortized capital cost rate of component r, 
including fixed charges and maintenance. A specific computer 
program was used for the numerical solution of the optimization 
problem consisting of the main program, double precision 
function, constraints subroutine and the optimization algorithm 
developed by Frangopoulos himself. In the statement of the 
functional optimization problem, the same objective function 
(Eq. (27)) was used. In the two methods, the same independent 
variables were taken into account as in the previous applications of 
conventional optimization method [36] and exergoeconomic 
optimization [37]. In the solution of the functional optimization 
problem, Lagrange multipliers were used. The simplified form of 
the Lagrangian was defined for the CGAM problem as [38] 


ArY,) + (M02 — A02¥o2) + A3¥3.0 + A5950 (28) 
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where 
6 

Py = Zp + ArrYrsr, r= 1,2,...,6 (29) 
T 


The required conditions 


VzL(x,y,à)=0, VyL(x,y,à)=0, V,L(x,y,à)=0 (30) 


were imposed on the Lagrangian to obtain the equations 


ə 6 
rÈ =0 (31 
r= 
Əlo2 
ho2 =a = C€ 32 
j ƏYo2 d ( 
_ Oly 
Àr Jy, (33 
Ayr = Ar (34 


In the solution of the nonlinear system presented for the CGAM 
problem, numerical methods were used and the minimization of 
Lagrangian (Eq. (28)) with respect to x by means of a nonlinear 
programming algorithm that was found to be more efficient than 
the analytic derivatives (by avoiding infeasible points) procedure 
[38]. The modular approach to the CGAM problem was given for 
demonstration purpose in the paper since the system presented 
through the problem is not very complex and detailed design 
characteristics are not necessary. The comparable optimization 
results for the nominal values of parameters in terms of three 
methods presented [38] are given in Table 6. 

Among three methodologies applied by Frangopoulos in his 
paper as part of sequential studies of the CGAM problem in the 
same issue, the direct use of nonlinear algorithm was the simplest 
one because it required the least effort in system analysis, but it did 
not include any information about the internal economy of the 
system (i.e. physical and economic relationships among the 
subcomponents of the CGAM system) [38]. The CGAM problem 
was solved by TFA successfully since the distribution of functions 
developed for subcomponents of the system established inter- 
relations between components or between system and the 
environment leading to a functional diagram of the system 
[22,23]. However the capabilities of TFA and IFA which was a 
further development of TFA are better revealed in more 
complicated situations than CGAM problem, when the solution 
of the complete optimization problem is required [22-25,38]. The 
modular approach has some significant advantages when detailed 
design of components is required. This methodology can be 
combined with direct nonlinear algorithms and TFA or IFA method; 
in particular decomposition of the optimization problem is 
applied. In Frangopoulos paper for CGAM problem, this was not 


Table 6 

Comparison of the nominal values of the decision variables and total cost of the 
plant for the CGAM problem through direct use of the nonlinear programming 
method, thermoeconomical functional approach and modular approach [38]. 


Decision variables Method 


Direct use of the TFA 
nonlinear algorithm [38] 


Modular approach 


PP, 8.59730 8.59770 8.59050 

Nsc 0.84641 0.84650 0.84653 

Nst 0.87886 0.87871 0.87878 

Ts (K) 912.77 913.14 912.93 

Ta (K) 1491.40 1491.97 1491.50 

F ($/year) 1.0426 x 107 1.0426 x 10’ 1.0426 x 10” 


the case, and therefore modular approach was applied to the 
problem only with demonstrative purposes. 


5.4. Engineering functional analysis approach (EFA) 


von Spakovsky constructed a thermoeconomic structure for the 
CGAM problem by applying engineering functional analysis 
method (EFA) [39]. Methodology in his paper consisted of two 
distinct parts: an economic model which represents the thermo- 
dynamic and economic behavior of the system and an algorithm to 
solve the model for optimum. The general form of objective 
function for CGAM problem in this paper was defined as [39] 


ee Soo S 8 28 8 3 
minimize T'ota = Ce + Capu + Coc + Cr + Curse + L puet + K 


(35) 
with respect to 


x = {Xi} = {Xc}, 


subject to system products, component inputs, junctions, 
branches, and physical limits respectively as 


i=1,2,...,I and c=C,APH,CC,T,HRSG (36) 


hpj=Y;j Y ;(Xpj) = 0, j=1,2...,P (37) 


he = Yj — Y;(xc Y) =0, j=P+1...,C and c 


= C, APH, CC, T, HRSG (38) 


2 
him =Yim — X Yjm=0, j=C+1...,L and m=1...,L (39) 
in=1 


2 
hen = >> Yenou-Yj=0, j=L+1...,B andn=1...,B (40) 


out=1 
Z(xX,y)>0, k=1...,K 41) 
where 
Ë = C (xe, Yc), c = C, APH, CC, T, HRSG 42) 
T na =T Suel(Xcc, Ycc) 43) 


Each of the general mathematical functions defined in Eqs. (35) 
through (43) were derived based on the equations given in the 
introductory paper of CGAM problem series [36,39]. The internal 
economy of the system was determined by EFA using Lagrange’s 
method of undetermined multipliers since with this approach, the 
real unit or marginal costs associated with developed thermo- 
economic structure and model in the paper were calculated. Thus 
Lagrangian was constructed from the objective function and 
equality constraints such that [39] 


. no -$ R 7 
Ltotal = DAC F L hiel = Vai = Y (Xp, )] T 5 Al; 
c=1 j=l 


j=P+1 
L 2 
-_ Y j(Xc, Y¢)] _ 5 Àj (v. = YY na} 
j=n+1,m=1 in=1 
B 2 
- > u( XO Ysp — r) (44) 
j=L+1,n=1 out=1 


The methods used in the optimization process in terms of EFA 
method were categorized in the paper as: modular, the Lagrangian 
and decomposition approaches. For decomposition approaches, 
two different sets of subgroups were designated namely Case 1 and 
Case 2 [39]. In Case 1, the first set of subgroups divided the system 
into two subgroups: the air compressor as one and the remaining 
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Table 7 
Optimum values of the objective function, component and fuel costs and the 
decision variable set for different methodologies of optimization through EFA [39]. 


Design Modular Lagrangian Decomposition Decomposition 
variables approach approach approach (Case 1) approach (Case 2) 
P2/P, 8.59822 8.59859 8.49996 8.49996 
se 0.84678 0.84675 0.84991 0.84991 
Nst 912.601 912.281 924.764 900.09 
T (K) 0.87874 0.87900 0.88000 0.87991 
T4 (K) 1491.566 1491.440 1528.970 1480.153 
Fora (S/S) 0.36204 0.36204 0.36850 0.36336 
Čat (S/S) 0.03653 0.03658 0.04740 0.03467 
T fuet ($/S) 0.32552 0.32546 0.32110 0.32869 
components as the other 
P2 
Xsub1 = Pi’ Nic and Xsub2 = {T3, T4, Ny} (45) 


In Case 2, the second set of subgroups divided the system into 
four groups: the air compressor as one, the air preheater, 
economizer and evaporator as another, the combustion chamber 
as third and the gas turbine as fourth. Its decision sets were 


P. 
Xsub1 = Pe me} Xsub2 = {T3}, Xsub3 = {T4}, Xsub4 


= {1} (46) 


A comparison of some of the major results obtained with each 
approach in terms of EFA method is given in Table 7. Optimum 
results for modular and Lagrangian were obtained as almost 
identical because these two methods coincided at the optimum 
(see Table 7) [39]. The decomposition approaches for both cases 
had some significant differences. This was expected since the 
decision variables sets defined for each subgroup were not disjoint. 
In Case 1, more realistic results were obtained compared with Case 
2 although four subgroups decomposition solution would be 
expected better than two subgroups solution. The analysis and 
optimization methodologies applied through this paper showed 
that EFA methodology is a very general one and therefore very 
useful when applied to the modeling, optimization and analysis of 
energy systems either time dependent or not [39]. 


5.5. Exergetic cost theory (ECT or TEC) approach 


Exergetic cost theory (ECT) [33] formulation of the optimization 
problem for CGAM system was based on the use of Lagrange 
multipliers [40]. The optimization problem, using the structural 
information provided for CGAM problem in the paper was 
formulated in a compact form as 


minimize T = ‘c.F + ‘uZ (47) 
subject to constraints 


KpP = F(I), Ps + (PF)F = P(II), 


k= f(t), r= fiT) @= fa(T) (49) 


where JT is the total cost of the plant product, Ce is a vector which 
contains the unit price of the resources entering the plant; F, and P, 
are the vectors containing the exergy of the entering resources and 
of the products leaving the plant respectively; u is a unit vector; Zis 
the vector containing the cost of the operating and maintenance of 
each component; Kp is a diagonal matrix containing the unit 
exergy consumption of each subsystem; P and F are vectors 
containing the exergy of the product and of exergy resources 
consumed (fuel) of each component; and (PF) is a matrix which 
defines the productive interconnection of the components [33,40]. 


P; = datum(III) (48) 


Table 8 
Optimization results for decision variables and costs in the 
CGAM problem by ECT methodology [40]. 


Decision variables Optimum values 


P2/P, 8.5234 

Msc 0.8468 

Nst 0.8786 

T3 (K) 914.28 

Ty (K) 1492.63 
Total cost rate ($/s) 0.36201 

Fuel cost rate ($/s) 0.32549 
Investment cost rate ($/s) 0.03652 
Cost of combustion chamber ($) 0.1469 x 10° 
Cost of air compressor ($) 0.1348 x 10” 
Cost of gas turbine ($) 0.1927 x 107 
Cost of air preheater ($) 0.8277 x 10° 
Cost of HRSG ($) 0.1202 x 107 


The optimization problem with constraints was transformed into 
one without constraints as 


minimize L(X,A) = M(X) + X- A;G;(X) (50) 


where /Tis the original objective function of the CGAM problem, Xisa 
vector containing the variables of the problem and A, is the Lagrange 
multiplier associated with the constraint G;. It was demonstrated in 
the paper that Lagrange multipliers associated with the constraints 
had a profound economic significance [40]. Optimization results of 
the CGAM problem by ECT is given in Table 8 [40]. 

The structural information obtained at the optimum through 
the ECT method is very useful when carrying out a thermo- 
economic analysis of the plant. Since CGAM problem is not a 
complex type of problem, the structural information developed 
through the analysis and optimization can be inferred from 
conventional thermodynamic analysis. In complex energy systems 
it would be very difficult [40]. 


5.6. Modified productive structure analysis (MOPSA) approach 


Kim et al. [45] proposed and applied a thermoeconomic method 
called modified productive structure analysis (MOPSA) to the CGAM 
system and emphasized how the cost structure of the cogeneration 
system was affected by the chosen level of aggregation which 
specifies the subsystems. Kwak et al. [60] investigated the cost 
structure of this predefined cogeneration system with the method 
proposed by Kim et al. [45] and he rewrote the exergy balance 
equations by reflecting the exergy losses due to heat transfer. Thus, 
assigning a unit exergy cost to every decomposed exergy stream in 
the cogeneration system, non-adiabatic exergetic cost-balance 
equations were written as [60] 


«CHE «BQ «BQ FY i iT 
Ex Cot io Ew = 5 Ex Ceo + (se, = 5 Eu) Cr 
input out put inlet outlet 
g , ; a Ġ& 
A (se — 5 é) Cp + To (xs = 5 Sj +) C; 
inlet outlet inlet outlet 
+Ż = Ė"Cw (51) 


The exergetic cost-balance equation given in Eq. (51) yields the 
productive structure of the thermal system [45,51] that was 
developed by Lozano and Valero [33]. Kwak et al. [60] described 
three levels of aggregation cases in their paper. In the first case (Case 
I) which was the simplest one, there was no internal parameter to 
predict the unit cost of products, and it yielded the lowest unit cost of 
electricity and the highest unit cost of steam. The second case (Case 
II) was the next simplest one. However, in this case there was only 
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Table 9 

The unit costs and the cost flow rates of the electricity and steam of the CGAM 
system by applying MOPSA approach in three different levels of aggregation cases 
[60]. 


Cost structure Case | Case Il Case Ill 
Unit cost of electricity ($/GJ) 8.12 8.58 8.46 
Cost flow rate ($/h) (879.96) (926.64) (913.68) 
Unit cost of steam ($/GJ) 11.62 10.55 10.83 
Cost flow rate ($/h) (533.24) (484.14) (496.99) 


one internal parameter, the unit cost of negentropy (i.e., negative 
entropy which is defined as the entropy that the system exports to 
maintain its own entropy level low), and with this level of 
aggregation case the highest unit cost of electricity and the lowest 
unit cost of steam were obtained. Case III represented the highest 
level of disaggregation. It had three internal parameters; the unit 
costs of thermal and mechanical exergy and negentropy. This 
yielded average results between Case I and Case II. In Table 9, 
comparative results of three aggregation levels for the CGAM 
problem are given. The production costs for the CGAM listed in 
Table 9 satisfies the following overall cost-balance equation for the 
system. Note that the unit cost of electricity and steam cannot be 
evaluated from this equation [45,51,60]. 


«CHE aP 5 3 i 
E, Cot+tE,Cu+ > 2 = WreeCw + IÉ" [Ceo (52) 


Table 11 


Table 10 
Comparison of exergy costing methodologies for the most disaggregated CGAM 
system in terms of the unit costs of electricity and steam. 


Methodology SPECO Structural MOPSA Exergetic 
[37] theory (STT) (Case III) cost theory 
[34] [60] (ECT) [33] 
Unit cost of electricity ($/GJ) 7.80 7.42 8.46 7.55 
Cost flow rate ($/h) (842.40) (801.36) (913.68) (815.18) 
Unit cost of steam ($/GJ) 10.45 10.97 10.83 10.64 
Cost flow rate ($/h) (479.60) (503.46) (496.99) (488.06) 


In Table 10, various exergy costing methodologies for the most 
disaggregated CGAM system are compared in terms of unit costs of 
electricity and steam. The unit costs of products of CGAM system 
obtained from exergetic cost theory (ECT) and SPECO approach 
[33,37] were close to the values of the simplest aggregation level 
(Case I). 

As it was confirmed in related works [33,34,37,51,60], the unit 
cost of products is dependent on the level of aggregation of the 
system when the irreversibility occurred in subsystems is given 
explicitly in the cost-balance equations. This situation suggests 
that the cost structure of the system is mostly affected by the 
irreversibility represented by the entropy production rate at each 
component, and that cost structure of the system is interrelated by 
the irreversibility occurred at each component [60]. 


Some major works on the thermoeconomic analysis and optimization of combined heat and power production systems. 


Author(s) 


Thermoeconomic methodology 


Type of system 


Tsatsaronis et al. [1] 
Tsatsaronis and Ho-Park [4] 
Frangopoulos [22] 
Frangopoulos [23] 
Frangopoulos [24,25] 
Frangopoulos [26,27] 

von Spakovsky and Curti [29] 
von Spakovsky and Evans [30] 
Tsatsaronis and Winhold [31] 
Lozano and Valero [33] 
Erlach et al. [34] 

Valero et al. [36] 

Tsatsaronis and Pisa [37] 
Frangopoulos [38] 

von Spakovsky [39] 

Valero et al. [40] 

Lazzaretto et al. [41] 

Verda [42] 

Verda and Borchiellini [43] 
Hua et al. [44] 

Kim et al. [45] 

Cerqueira and Nebra [46] 
Manolas et al. [47] 
Lazzaretto and Tsatsaronis [48] 
Kwon et al. [51] 

Lazzaretto and Tsatsaronis [52] 
Cziesla and Tsatsaronis [54] 
Kwak et al. [57] 

Kwak et al. [60] 

Valdes et al. [59] 

Silveira and Tuna [62,63] 
Vieira et al. [65,75] 

Zhang et al. [67] 

Erdil [68] 

Mazur [69] 

Bonnet et al. [70] 

Colpan and Yesin [73] 

Li et al. [76] 

Ozgener et al. [77] 
Zaleta-Aquilar et al. [78] 
Unver and Kilic [79] 
Abusoglu and Kanoglu [80,81] 


LIFO 

SPECO/AVCO 

TFA 

TFA 

IFA 

TFA and optimization 
EFA 

EFA 

EEA 

TEC 

STT 

Conventional solution 
SPECO/AVCO 

TFA 

EFA 

TEC 
Thermoeconomic diagnosis (TD) 
TD 

STT and TD 

EEA 

MOPSA 
SPECO/TFA/TEC/TECD 
TFA and genetic algorithm (GA) 
EEA 

SPECO/MOPSA 
SPECO/EEA 

IEA— fuzzy logic 

EEA 

MOPSA 

TFA and GA 

EEA and optimization 
STT 

STT 

EEA and optimization 
EEA- fuzzy logic 

EEA 

SPECO 

TFA and optimization 
SPECO 

EEA 

TEA 

SPECO 


Simple gas turbine cogeneration 
Simple gas turbine cogeneration 
Complex thermal system 

Steam cycle cogeneration 

Steam cycle cogeneration 

Steam cycle cogeneration 
Cogeneration/heat pump facility 
Combined cycle cogeneration 
Combined cycle cogeneration 
Simple thermal power plant 
Combined cycle cogeneration 
CGAM problem 

CGAM problem 

CGAM problem 

CGAM problem 

CGAM problem 

TADEUS problem 

TADEUS problem 

TADEUS problem 

Combined cycle cogeneration 
Gas turbine cogeneration system 
Simple gas turbine cogeneration 
Industrial cogeneration system 
Combined cycle cogeneration 
Gas turbine cogeneration system 
Combined cycle cogeneration 
Thermal power plant 

Thermal power plant 

CGAM problem 

Combined cycle gas turbine plant 
Combined cycle cogeneration 
CGAM problem 

Coal fired thermal power plant 
Combined cycle cogeneration 
Combined cycle cogeneration 
Micro-cogeneration 

Combined cycle cogeneration 
Dist. multi-generation systems 
Geothermal district heating 
Steam turbine 

Combined cycle cogeneration 
Diesel engine cogeneration 
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6. Exergoeconomic analysis and optimization of combined 
heat and power production (CHPP) 


Combined heat and power production (or cogeneration) is the 
simultaneous production of power and usable heat from the same 
fuel or energy source [77,79-81,84]. Cogeneration typically 
improves the overall efficiency of energy use, leading to reduced 
fuel use and pollutant emissions, and is an example of distributed 
generation [84,85]. Depending on the purpose, exergoeconomic 
studies on combined heat and power production in literature can be 
grouped in two categories: (1) if the purpose is the allocation and 
determination of the production costs and/or cost of losses of the 
plant, the study is focused on the monetary flow rate accounting 
through components in the plant [4,31,48,57,66,67,70,73,77- 
81,84]. (2) If the purpose is the optimization of a power system, 
the study is focused on the selection of the best conditions for 
operating the system [23-29,36-45,47,50,51,53-55,59-63,64, 
68,75]. Thermoeconomic concepts and analysis and optimization 
of combined heat and power production have been studied in a 
satisfactory amount of works. Among these publications, some 
specific and remarkable ones are given together with the thermo- 
economic method used and type of system studied (Table 11). In the 
published papers presented in Table 11, several thermoeconomic 
methods were applied to obtain the unit costs of the products for 
combined heat and power systems. Application of different methods 
to a given system normally yields different values of product cost 
because the assumptions made in the formulations of cost-balance 
equations are different in each method. Besides, the cost structure of 
combined heat and power systems is affected by the chosen level of 
aggregation that specifies the subsystems. The level of aggregation 
depends on the number of parameters performed in the determina- 
tion of the unit costs of products [2]. Note that if the irreversibilities 
in system components are not considered the unit cost of the 
products does not depend on the chosen level of aggregation [60]. 

A major fraction of published works listed in Table 11 were 
performed using algebraic methods [1,4,31,33,37,45,46,51,52,54, 
60,73,77,80,81]. In these studies, SPECO and MOPSA approaches 
were applied to combined heat and power systems in order to 
reduce the subjectivity of fuel and product definitions and cost 
partitioning. This is done by the help of an iterative exergoeco- 
nomic performance improvement procedure based on exergoe- 
conomic variables such as relative cost difference, exergoeconomic 
factor and exergetic efficiency as introduced by Tsatsaronis 
[1,2,4,31,33,37]. 

The remaining studies were performed using calculus methods 
[22-27,29-34,38,39,44,48,59,70,76,78] mainly developed by Fran- 
gopoulos [22-27] and von Spakovsky [29,30,39]. As mentioned 
before, these methods are based on Lagrange multipliers and they 
vary from iteration to iteration when thermoeconomic isolation of 
component is not achieved. This was the main difficulty of the 
method in the application. After the introduction of functional 
diagram concept and by assigning only one function and one 
product to each component by Thermoeconomic Functional 
Analysis (TFA) method [22-27], the complex nature of the calculus 
methods were overcome. The Lagrange multipliers show relatively 
smaller variations in the optimization process with the application 
of second law in the TFA method [60]. 

In order to reduce the complexity in the optimization of 
complex thermal systems, several decomposition methods based 
on the second law were proposed [44,54,62,63,69]. Combined heat 
and power systems have some large scale problems due to their 
complicated nonlinear characteristics. Thermodynamic restric- 
tions make the optimization problem difficult to be solved. Many 
researchers agree that focusing on the dominant decision variables 
is more important than manipulating all variables simultaneously 
in the optimization case [60]. 


7. Conclusions 


Thermoeconomics assesses the cost of consumed resources, 
money and system irreversibilities in terms of the overall 
production processes. Methodologies presented through thermo- 
economics may help to point out how resources may be used 
effectively for sustainable development. Assessing the cost of the 
flow streams and processes in a complex system helps to 
understand the process of cost formation from the input resources 
to the final products. These analyses can solve problems related to 
complex energy systems such as combined heat and power 
production (CHPP) that could not be solved by using conventional 
energy analysis. 
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